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ABSTRACT 

We present the results of the analysis of a sample of 17 low-luminosity 
(-^X ;$ 10 42 erg/s), radio loud AGNs in massive galaxies. The sample is extracted 
from the SDSS database and it spans uniformly a wide range in optical [OIII] emis- 
sion line and radio luminosity, but within a narrow redshift range (0.05 < z < 0.11) 
and a narrow super massive black hole mass range (~ 10 8 M Q ). For these sources we 
measured core X-ray emission with the Chandra X-ray telescope and radio emission 
with the VLA. Our main goal is to establish which emission component, if any, can be 
regarded as the most reliable accretion/jet-power estimator at these regimes. In order 
to do so, we studied the correlation between emission line properties, radio luminosity, 
radio spectral slopes and X-ray luminosity, as well as more complex multi-variate rela- 
tions involving black hole mass, such as the fundamental plane of black hole activity. 
We find that 15 out of 17 sources of our sample can be classified as Low-Excitation 
Galaxies (LEG) and their observed properties suggest X-ray and radio emission to 
originate from the jet basis. We also find that X-ray emission does not appear to be 
affected by nuclear obscuration and can be used as a reliable jet-power estimator. More 
generally, X-ray, radio and optical emission appear to be related, although no tight 
correlation is found. In accordance with a number of recent studies of this class of ob- 
jects these findings may be explained by a lack of cold (molecular) gaseous structures 
in the innermost region of these massive galaxies. 
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1 INTRODUCTION 

In this work we focus on three accretion/jet-power estima- 
tors among the different signatures of galactic nuclear ac- 
tivity arising at different wavelengths: optical line emission, 
nuclear X-rays and nuclear non-thermal radio emission. In 
the standard "unified model" for AGNs (|Urrv fc Padovaml 
1 19951 ). optical narrow emission lines come from gas located 
several hundreds pc away from the central engine. These 
lines are excited by ionizing radiation produced in the in- 
nermost accretion flow and escaping along the polar axis of 
the obscuring torus that surrounds the black hole. Since this 
ionized gas is so distant from the central engine, the obscur- 
ing torus does not affect greatly its flux, thus narrow optical 
emission lines suffer only moderate amounts of dust obscu- 
ration due to the interstellar medium. This suggests that 
[OIII] emission line luminosity can be a good estimator for 
the AGN accretion power. On the contrary, X-ray emission 



arises directly from the hot corona surrounding the accre- 
tion disk or from the base of a relativistic jet (for radio-loud 
objects). As such, it represents a more faithful estimator of 
the accretion power, although it could be heavily obscured 
by high column density through the dusty torus in objects 
where we look at the AGN from a close to edge-on sight-line. 
Finally, powerful non-thermal radio emission is the observa- 
tional signature of the presence of a jet whose relativistic 
particles emit synchrotron radiation going through strong 
magnetic fields. Radio emitting jets are ubiqu itous, particu- 
larly at low intrinsic powers (see e.g. lHoll200Sl , and references 
therein), and it has been postulated that it also could be 
used to estimate the AGN kinetic and tot al power output 
jFalcke fc Biermannll 19951 : iBest et~al]|2006h. 



As emerged in recent years (see lHardcastle et alj|2009l . 
and references therein), radio loud objects can be divided 
into two main populations: 

high excitation galaxies (HEG): quasars and high-power 
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narrow line radio galaxies that are likely to be the same 
population, seen at different orientation; 

low excitation galaxies (LEG): active radio galaxies of low 
intrinsic power characterized by emission line spectra of low 
excitation, that probably form the parent population of the 
mostly lineless BL-Lac objects. 



Recent studies ijChiabe rgc ct al. 



Whvsong fc Antonucci 12004 lEvans et all 



2002; 



2006; 



Hardcastle et al.1 I2006I . 120091 ) have proposed that LEG 



lack any of the conventional apparatus of the AGNs 
(i.e. radiatively efficient accretion disk, X-ray corona and 
obscuring torus) and that their radio and X-ray nuclear 
emission can be explained as a result of the properties of 
the small-scale jet. 

Significant progress has also been achieved in un- 
derstanding how the radio properties of AGN relate to 
other paramet e rs such as bla c k hol e mass and accr etion 
rate llHcll2002h. iMerloni et all (|2003l ) (hereafter |M03| ) and 
iFalcke et al.l (|2004h considered heterogeneous black hole 
samples spanning a broad range in black hole mass and 
X-ray luminosity, and found a strong correlation between 
the radio luminosity Lr from the unresolved radio cores of 
AGNs, the black hole mass M and the 2-10 keV X-ray lumi- 
nosity Lx of the form (|M03l ): 



log Lr = 7.44 + 0.6 log L x + 0.78 log (M/M ) 



(1) 



These scalings with mass and luminosity can be well ex- 
plained by the theory of synchrotron emitting compact 
relativistic jets and radia t ively inefficient accretio n flows 
l|Blandford fc Konigll [l979l : iHeinz fc Sunvaevll2003h . These 
works explicitly demonstrated that the black hole mass 
should be considered as a fundamental parameter in the 
determination of the observed relations between emission 
components in various bands. 

In this paper we want to push forward these investiga- 
tions on the relations between different accretion and power 
output estimators in super massive black hole (SMBH). In 
particular, we will test whether the loose empirical [OIII] 
- X-ray correlation found in type 1 AGN, as well as the 
fundamental plane of black hole activity, hold in an unbi- 
ased sample of type 2 AGN in the local univ erse. Compared 
to pr evious wo rk on this subject (see e.g. He ckman et al.l 
2005, hereafter H05), this study is unique because we select 
objects in a very narrow range of black hole masses, thus 
disentangling the effects of mass and accretion rate on the 
observed [OIII] - X-ray - radio correlations. 

The paper is organised as follows: in the next section 
we describe our sample, while in Section [3] we lay out the 
most important calibration steps used for the X-ray, radio 
and optical data reduction. Then in Section [4] we examine 
the relation between different power estimators. Finally in 
Section [5] we discuss their implications. 



2 SAMPLE SELECTION 

In order to answer some of the important questions posed 
in the previous section, we have performed a snapshot X- 
ray/radio survey of a new, well defined sample of low- 
luminosity AGN. The sample used for our analysis was 
originally extract ed from 2712 radio luminous type 2 AGN 
\Best et al.l |2005( ) contained in the 212000 galaxies of the 



second data release of the SDSS. From this sample we se- 
lected 17 objects with the following characteristics: 

• The black hole masses have been deri ved from veloc- 
ity dispersion measuremen ts of each galaxy l|Tremaine et al.l 
2002; Heck man et al.ll2004f ). We require black hole masses to 
be 8 < log A/bh < 8.5 (the range is comparable to the un- 
certainty in the estimate of A/bh). This particular window 
was chosen to maximize the number of objects with z < 0.11 
that are radio loud. 

• A second cut was made on the redshift: 0.05 < z < 0.11. 
The SDSS spectra are obtained using 3 arcsec diameter fi- 
bres and at larger distances, the AGN emission is signifi- 
cantly contaminated by starlight from the host galaxy. 

• The objects were chosen to span a wide range in dust- 
corrected [OIII] luminosity (i[oni] ~ lO 5,5 Z/0 up to to 
Aoiii] = 10 9,2 £©), and the entire observed range in radio 
luminosity, 10 38 5 erg s _1 < Lr < 10 413 erg s" 1 . 

• In order to not be biased by cluster environment and 
to facilitate the X-ray data reduction, we chose objects that 
were at least 7 arcminutes away from the centre of galaxy 
clusters. 

142 AGN satisfy all these conditions. Because of their in- 
trinsic weakness, a full multi-wavelength survey of the entire 
sample was deemed too time-consuming. Instead, we chose 
to focus on a randomly selected sub-sample of 17 objects 
to cover the L[oiii] ~~ Lr parameter space widely. The fi- 
nal choice of the sample size was made by trading off the 
total observing time requested for maximizing the chances 
of source detection in the X-rays and ensuring that enough 
sources are included for statistical tests. 

Properties of the galaxies in the selected sample are 
listed in Table Q] while in Figure [T] we plot the sample on 
the [OIII] - radio plane using line luminosities from SDSS 
data release 7. At the time of the sample selection the last 
SDSS release available was the fourth. Subsequent analysis 
and better calibration of the SDSS spectroscopic data im- 
ply that some [OIII] luminosity values differ slightly between 
the two releases due to different data reduction procedures 
and this is the reason why two objects (9 and 17) are quite 
near in the [OIII] - radio plane while, as stated above, we 
selected them to cover that space as uniformly as possible. 
Finally, we stress that in the sample selection process we did 
not exclude any source on the basis of the ionization state of 
the emission line complex, so that no bias for or against the 
HEG/LEG nature of the sou rces had been introdu ced. After 
our observations were taken, IPlotkin et al.l (J2008) , using ra- 
dio and optical data, have classified one of them as a BL-Lac 
with a strong boosting in the radio and X-ray emission. This 
is source number 4 in Table [1] and is a clear outlier in most 
of the diagrams present throughout the paper. We decided 
to leave it in the sample for completeness of information. 

Throughout this paper luminosity distances are cal- 
culated adopting Ho = 70 km s _1 Mpc _1 , Q. m = 0.3 and 
n A = 0.7. 



3 DATA ANALYSIS 

Our sample of 17 objects has been observed with Chandra 
X-ray telescope (see Section 13. ip and with the VLA (see 
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Table 1. The sample 



Object ID RA DEC z logL[ OIII] logL x log^R. «h logM BH 

J2000 J2000 2-10 keV 1.4 GHz 4.8 GHz 8.4 GHz 
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NOTE: Col. (1): Name of the object. Col. (2): Identification number. Col. (3): Redshift. Col. (4): [OIII] luminosity log(L[ OIII j [erg/s]) 
corrected for dust absorption. Col. (5): X luminosity log(Lx[erg/s]) (corrected for absorption). Col. (6): Radio luminosity at 1.4 GHz 
log(Li.4 [erg/s]). Col. (7): Radio luminosity at 4.8 GHz log(L4.g[erg/s]) (mean 5 <r error: 0.023 dex). Col. (8): Radio luminosity at 8.4 
GHz log(Lg.4[erg/s]) (mean 5 a error: 0.024 dex). Col. (9): Radio spectral index, defined with F oc v~ a . Col (10): SMBH mass 

log(M[MQ]). Object signed with a f is resolved in the radio maps. 
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Figure 1. Sources in our sample plotted in a radio (1.4 GHz) - 
optical ([OIII]) luminosity plane. Sources are selected in order to 
cover wide range of luminosity both in the radio and the optical. 
Points are coloured according to their spectral index in the radio 
band (red: steep-spectrum — blue: flat-spectrum, see section 13.21 
for details). Square markers are for LEG steep-spectrum objects, 
round markers are for LEG flat-spectrum objects while the tri- 
angle marker is the HEG and the star marker is the BL-Lac. For 
the classification criteria see section [3.21 and 13.31 



Section [3.2[) . Information about their [OIII] luminosities is 
taken from SDSS and is analysed as explained in Section [3T3T 
All errors due to inaccurate evaluation of distances 
are considered negligible, while we account for an erro r 
of 0.25 dex in the black hole mass jTremaine et alj|2002r i. 
Throughout the paper linear fits have been performed us- 
ing the Buckley - James method l|Bucklev fc Jamesl Il979l : 
llsobe et alj Il986l ) to account for upper limits. With this 
regression algorithm, the most probable values of the up- 
per limits are evaluated through the Kaplan- Meier non- 
parametric distribution (jKaplan fe Meierl Il958l ) and then 
used as observed values for a linear fit that is performed us- 
ing the orthogonal distance regression to take into account 
errors both in the dependent and in the independent vari- 
able. Correlation tests we re performed usin g a generalized 
Kendall correlation test jlsobe et al.l [T986) that accounts 
also for upper limits an d, if needed, for partial correlation 
l|Akritas fc SieberrJ[l996l ). 



3.1 X-ray 

All of the AGN were observed on the S3 chip of the Ad - 
vanced CCD Imaging Spectrometer (ACIS: rCarmird 120031 ) 
aboard the Chandra X-ray Observatory. Data were taken 
in timed exposure mode with the standard frame time at 
the default location on the S3 chip and telemetered to the 
ground in very faint mode. Data reduction and point source 
extractions were completed using CIA Q versions 4 . 2, an d 
ACIS Extract version 2010-02-26 (AE; iBroos et alj|20ich . 
respectively. During extraction, the source position was ad- 
justed based on the mean position of the extracted counts 
to more accurately calculate the point source photometry 
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(except for the two sources with upper limits). The posi- 
tion refinement was always well within Chandra's absolute 
position uncertainty. 

We used Sherpa version 4.2 to jointly fit the unbinned 
source and background spectrum of each source dFreeman 

musing the C-statistic, which is similar to the ICash 
statistic but with an approximate goodness-of-fit 
measure, and the moncafQ minimization algorithm. We used 
this statistical method for the fitting to avoid losing the little 
spectral information that we have for most of our spectra, 
which frequently have very few counts fe.g-. lNousek fc Shud 
Il989h . We fit each source with an absorbed power-law model 
(xszphabs x xspowerlaw). In all cases, the background ex- 
hibited the signs of hard particle background and was fit by 
two power laws and thermal component, which achieved a 
good fit. These ad hoc models are only used to constrain the 
background component of the spectra. 

Since degeneracies in the fit parameters will frequently 
arise for very faint sources, we followed a specialized scheme 
for these sources based on the number of counts in the source 
extraction region. If we extracted less than 5 counts (0.5-8 
keV) for a source, the power-law index and the hydrogen 
column density for the fit were frozen to 1.7 and the Galactic 
value, respectively. If we extracted more than 5 but less than 
26 counts (0.5-8 keV), we froze only the power-law index to 
1.7 and let the hydrogen column density float, although, in 
this case, it was always poorly constrained. For all other 
sources with more than 26 counts (0.5-8 keV), we allowed 
all fit parameters to float. 

As can be seen in Table [2j the AGN span a wide range 
of brightness. We put upper limits on the luminosities of 
two of the sources with only two counts each (identified 
with number 2 and 5 in Table [T}. On the other hand, 
there were also two sources (IDs 4 and 7 in Table [1} that 
had non-negligible pileup fraction^ (39.3% and 4.6% of the 
counts, respectively) and were therefore fit with the stan- 
dard pileup model. Only two sources (number 8 and, to a 
lesser amount, 13) have significant additional intrinsic ab- 
sorption above that expe cted for the galactic foreground 
l|Dickev fc Lockm an 1990). However, the determination of 
the column density for all of our sources is highly uncertain. 

The unabsorbed X-ray luminosities listed in Table 1 
were calculated from each of the model fits and the uncer- 
tainties were calculated as follows. The unabsorbed luminos- 
ity is evaluated at each point in a three-dimensional grid 100 
points on a side, each dimension corresponding to the num- 
ber of undetermined parameters in the source model (the 
absorption, power-law photon index, and power-law normal- 
ization). Then the minimum and maximum luminosity is 
selected within the confidence interval for the appropriate 
change in statistic value (AC = 3.53 for o ne sigma and 
three interesting parameters; e.g., lAvni|[l976l ). The photon 
index is constrained to be 0.7-2.7 and the column density 
must be at least the galactic minimum in the direction of 
the source. In many cases, the uncertainties in the lumi- 
nosity estimated with this method are very different from 
those estimated f rom the counts using Bayesian methods 
l|Kraft et al.lll99ll ). 



Finally, we note that, since the Chandra observations 
were carried out over the span of one year (late 2006-early 
2008), variability between the time of X-ray observations 
and radio observations could affect some of our results. 



3.2 Radio 

Radio data were taken in 2007 using the VLA in A con- 
figuration with a total observation time of 10 hours. The 
calibration procedure has been performed using CASA soft- 
ware (version 3.0.1) developed at NRAC0- At the time of the 
observations, 12 out of 26 antennas were already equipped 
with EVLA receivers. The presence of two different types of 
receivers on the antennas caused different response from the 
correlator whether a mixed baseline (VLA-EVLA) or a non- 
mixed one (VLA- VLA and EVLA-EVLA) were processed. 
To solve the problem, a baseline-based calibration was per- 
formed using the CASA task blcal. Finally, some cycles of 
self-calibration have been done on all sources. 

We observed our sample in two different bands (4.8 
and 8.4 GHz) and from these flux values radio spectral in- 
dexes have been extracted. Sources with a spectral index 
(F oc v~ a ) higher than a — 0.3 are classified as steep- 
spectrum, while other sources are classified as flat-spectrum. 
Values from the FIRST survey (1.4 GHz) were not taken into 
account during the spectral index extraction process because 
source fluxes could have changed since FIRST observations. 
From the spectral analysis we find that ten out of seventeen 
sources of our sample have a flat spectrum, while in one case 
(source 13), the spectrum is inverted. 

To check the radio calibration procedure we computed 
the spectral index average values using the three different 
bands (1.4, 4.8 and 8.4 GHz) two by two. If the main driver 
for changes in the spectral shape between any two bands is 
time variability, then some of the sources would have an in- 
creased spectral index while some other a lowered one, with 
an average value of zero. This means that the three averaged 
spectral index values should be similar. We find an averaged 
spectral index 04.8-8.4 = —0.29 ± 0.12 between 4.8 and 8.4 
GHz, 0:1.4-4,8 = -0.24 ±0.10 between 1.4 and 4.8 GHz and 
ax. 4-8.4 = -0.25 ±0.09 between 1.4 and 8.4 GHz. They are 
all compatible, enforcing the idea that individual differences 
between spectral indices including the low-frequency FIRST 
observations are mostly due to time variability. 

The sample was selected to be unresolved at the FIRST 
resolution (5") and only one source (number 17) shows clear 
evidence of a resolved structure in our new observations, 
that is a 1 arcsec long, jet-like feature in the west-north- 
west direction. Finally radio fluxes were corrected to be all 
at the same rest-frame frequency using the derived spectral 
indexes. Errors on radio fluxes have been set to 5 times the 
background RMS obtained analysing final maps. 



3.3 Optical 

We have re-analysed the available SDSS spectra for the 
sources in our sample in order to assess the physical na- 
ture of the emission lines. Optical spectroscopic informa- 
tion plays a major role in understanding emission proper- 



1 http:/ /cxc. harvard. edu/sherpa4.2/ahelp/montecarlo.py.html 

2 http:/ /cxc. harvard.edu/ciao/ahelp/acis_pileup. html 



3 http://casa.nrao.edu/ 
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Table 2. X-ray data analysis values 



ID 


ID 


Obs. 


Exposure 


Source 


Exp. bkg 


Photon Index 


N H [10 22 cm- 2 ] 


Probability 




Chandra 


date 


time [s] 


counts 
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(1) 


(2) 
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(6) 


(7) 


1 


8244 


2008-01-16 


11629 


12 
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0.00 


7 


8250 


2007-06-29 


6880 


260 
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2007-03-07 
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NOTE: Col. (1): Identification number. Col. (2): Identification number used in Chandra observations. Col. (3): number of source counts 
(0.5 - 8 keV) in the 90% PSF. Col. (4): expected background counts in source region (0.5 - 8 keV). Col. (5): is the best fits for the 
photon index and relative upper/lower limit increments. When the value is frozen to 1.7 (see text for details) the upper and lower 
limits are not calculable and therefore omitted. Col. (6): is the column densities Nh and relative upper/lower limit increments, f 
indicates that Njj is frozen to the galactic minimum. { indicates that the best fit happens to be at the galactic minimum. Col. (7): 

"prob_no_source" value, defined in section 5.10.3 of the AE user manual 
http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide/). Sherpa's "conf" routine has been used to calculate the uncertainties in 
the photon index and column density (col. 5 and 6). Sherpa's "conf" method gives similar uncertainties to the grid method used in the 
calculation of the uncertainty in the luminosity, except that Sherpa's "conf" routine tends to underestimate the uncertainty relative to 
the grid method (e.g., on the order of a few tens of percent for the photon index). 
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Figure 2. Diagnostic diagrams for our 17 sources. Black solid lines ijKewlev et al.|| 2006T) separat e star forming galaxies, on the bottom 
left of each diagram, from AGNs. The regio n between dashed an d solid lines llKewlev et al.ll2006f ) identify composite galaxies with both 
star forming and AGN activity. Red lines llKewlev et alj|2006h separate Seyfert galaxies (above) from LINERs (below). Green lines 
(Buttiglionc ct al. 200§) separate HEG (top) from LEG (bottom) sources. Average error bars arc shown in the top left corner of each 
plot. 



ties of central engine of radio galax ies. A firs t class ifica- 
tion using spectral li nes was made bv lHeckmanl £l980) and 
iBaldwin et al.l (|l98ll ). They proposed to use optical lines ra- 
tios to distinguish between HII regions ionized by star form- 
ing activities and regions ioni zed by AGNs. Some years later 
IVeilleux fc Osterbrockl (|1987m revised these definitions using 
only ratios of lines with a small separation in wavelength, in 



order to reduce the reddening problem and the impact of flux 
calibration errors. They introduced these lines combinations: 
[OIII]A5007/H/3, [NII]A6583/Hq, [SII]AA6716, 6731/Ha and 
[OI]A6364/Ha. The separat ion between AGN a nd HII re- 
gio ns was also ca l ibrate d bv lKewlev et al. (2001) and again 
bv lKewlev et al l l|2006h using ~ 85 000 emission line galax- 
ies from the SDSS catalogue. Seyfert and LINERs occupy 
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in these diagnostic diagrams different piaces and the au- 
thors suggested that the observed dichotomy correspond to 
two AGN subpopuiations associated with different accretion 
modes. 

Attempts to appiy thes e diagnostic d iagram s to radio 
ioud gaiaxies were made by lLaing et al.l l)l994h on a sub- 
sample of 3CR radio- galaxies. Foiiowing the suggestion by 
Iffine fc Longairl (|l979f ) that FRII sources can be divided into 
two subciasses, they separated their sample into high excita- 
tion galaxies (HEG, defined as galaxies with [OIII] /Ha > 0.2 
and equivalent width of [OIII] > 3 A) and low excitation 
galaxies (LEG). Following this idea, we have divided our 
radio lou d sample in HEG and L EG using the criteria de- 
scribed in lButtiglione et alj {2009.) 5 where they take advan- 
tage from the combination of the optical diagnostic planes. 
They define LEG all the sources with log([OIII]/H^) - 
[log([NII]/H Q ) + log([SII]/H Q ) + log([OI]/H a )] /3 < 0.95. 

We calculated lines luminosities using the spectra ex- 
tract ed from the SDSS data release 7 (|Abazaiian et al] 
2009). For our purposes we are interested both in the [0111] 
observed line fluxes and in the dust corrected line fluxes. 
When dust attenuation needs to be taken into account we 
measure Ha and H/3 line intensities and use the Balmer 
decrement method . We assume an intrinsic Ha/ H/3 of 2.87 
l|Osterbrocklll989r ). The attenuation correction we use is a 
double power law l|Charlot fc Fallll2000r ) of the form: 

s^x^r 3 ^^ ' 7 (2) 

where r„ is the total effective optical depth in the V band. 
With this double power law we account for the attenuation 
due the presence of discrete clouds random distributed (first 
term) and the attenuation due to the ISM (second term), so 
/i is the fraction of total r„ caused by ambient interstellar 
medium. We set /i = 0.3 based on observed relations between 
UV continuum slope and the Ha to Hj3 line ratios. 

We have generated emission line diagnostic (BPT, af- 
ter [Baldwin~^FaO[l98l]) diagrams to compare emission lines 
ratios ([OIII]A5007/H/3, [OI]A6003/Ha, [SII]A6717, 7631/Hq 
and [NII]A6584/Hq) of our sample (see Figure 0). 

Following th e prev iously explained classification from 
iButtiglione et al l (2009) we can classify the majority of our 
sources as LEG and be confident only for source 8 to be 
classified as an HEG. Some sources (11 and 12) fall in an in- 
termediate position since they do not satisfy requested line 
luminosity ratios in all diagnostic diagrams and they are 
then classified as LEG. We must notice that with this divi- 
sion all sources optically classified as LINERs are now clas- 
sified as LEGs and source number 8, optically classified as a 
Seyfert, is an HEGs. The difference in the classifications is 
mainly for borderline objects (in our sample: 11, 12 and 16) 
that are conventionally classified as Seyfert galaxies but in 
IButtiglione et al l (2009) classification are instead LEGs (see 
Figure HI). Contrary to what was found bv lButtiglione et al.l 
(2009) we have some strong [OIII] emitters (£[oiii] > 10 41 ' 5 ) 
that do not have log [OIII]/H/3 ~ 1, but lower values, e.g. 
source 12 and 14. Finally, object number 4 is an outlier in 
many plots. The rea son is that this source, as classified by 
iPlotkin et all (2008) using radio and optical data, is a BL- 
Lac with a strong boosting in the radio and X-ray emission. 



Table 3. Partial correlation coefficients 



Indep. 


Dep. 


Test 


Prob. of non 


correlation 


variable 


variable 


variable 


observed 


corrected 






All sampl 


e 




X-ray 


Radio 


[OIII] 


0.93 


0.98 


Radio 


[OIII] 


X-ray 


0.29 


0.06 




X-ray 


Radio 


0.01 


0.09 






LEG 






X-ray 


Radio 


[OIII] 


0.60 


n a a 

U.4o 


Radio 


[OIII] 


X-ray 


0.12 


<0.01 


[OIII] 


X-ray 


Radio 


0.23 


0.44 




Radio steep-spectrum sources 




X-ray 


Radio 


[OIII] 


0.38 


0.82 


Radio 


[OIII] 


X-ray 


0.22 


0.02 


[OIII] 


X-ray 


Radio 


<0.01 


0.09 




Radio fiat-spectrum sources 




X-ray 


Radio 


[OIII] 


0.45 


0.48 


Radio 


[OIII] 


X-ray 


0.87 


0.88 


[OIII] 


X-ray 


Radio 


0.05 


0.05 



NOTE: In boldface those values where the significance is > 95%. 
Object number 4 is not considered. 

4 COMPARISON AMONG [OIII], X-RAY AND 
RADIO PROPERTIES 

In this section we analyse the correlations and properties 
that connect [OIII], X-ray and radio emissions from the 
sample of sources described in Section [2] As a prelimi- 
nary step, we have performed a partial correlation analy- 
sis on the sample (and various sub-samples), both for the 
observed luminosities and for those corrected for absorp- 
tion/extin ction, as discussed in th e previous section. We 
adopt the lAkritas fc SieberiJ (l996) method to control for 
censored data. Table shows a summary of such a study, 
where the presence of a correlation between two variables is 
tested accounting for the common dependence on the third 
(test variable). In the next subsections relations between 
different accretion power estimators are evaluated and dis- 
cussed in detail: in section |4~T1 we present the analysis of the 
relationship between [OIII] and X-ray emission in our sam- 
ple, while in section 14.21 we exploit the knowledge of radio 
spectra to provide an interpretation of the X-ray and [OIII] 
information. Finally, in section 14.31 we discuss the X-ray - 
radio correlation (or lack of it) in light of previous studies 
on the fundamental plane of active black holes. 

4.1 Accretion estimators: the relationship 
between [OIII] and X-ray emission 

In |H05| the authors emphasized the strong difference that 
arises if one considers AGN samples of Seyfert galaxies se- 
lected by X-ray luminosity or using the [OIII] line luminos- 
ity. For a sample of hard X-ray selected AGNs (flux limit 
in the 3-20 keV band of 2.5 x 10" 11 erg s _1 cm -2 ), they 
find a mean log(Lx/ L[om]) of 2.15 dex and a = 0.51 dex. 
Type 1 and Type 2 sources are separated in terms of emis- 
sion intensities, but do not present significant difference in 
the mean ratio between X-ray and [OIII] luminosities. In 
the left panel of Figure 131 the lH05l X-ray selected dataset 
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Figure 3. Left panel: our sample o n a X -ray [2-10 keV] - optical ([OIII]) plane (uncorrected for absorption). Symbols are as in Figure[T] 
Black dots are Type 1 objects from lH05l hard X-ray selected sample, while circles are Type 2. Red and blue solid lines are the mean ratio 
between logLx an d logL[oin] calculated using respectively steep and fiat sources (coloured patches are at one sigma). Da shed lines are 
linear regressions (object number 4 is excluded). In the right panel the comparison is with the [OIII] selected sample from lHoR 



Table 4. log(Lx/£[oill]) f° r uncorrected data 



Sample 


log(L x /L[oni]) 


a 




[dex] 


[dex] 


H05 - X-ray selected 


2.15 


0.51 


H05 - [OIII] selected (Type 1) 


1.59 


0.48 


H05 - [OIII] selected (Type 2) 


0.57 


1.06 


This work 


1.25 


0.57 



is shown together with our sample, without applying any 
absorption/extinction corrections, for ease of comparison 
(absorption-corrected correlations are discussed in the fol- 
lowing sections) . Sources in our sample have a systematically 
lower X-ray to [OIII] ratio (see Table this is in agreement 
with t he di fferent selection criteria used for the two samples. 
In the Sol X-ray selected sample a higher X-ray emission is 
indeed expected. 

In the right panel of Figure [3] we compare instead our 
dataset with the [OIII]-selected sample of lH05l ([OIIIJA5007 
fluxes greater than 2.5 x 10 _13 erg cm -2 s _1 , corresponding 
to L [0 iii] > 10 40 erg s" 1 ). As noted by HoE Type 1 AGNs 
are detected at the expected rate, while around 2/3 of the 
Type 2 AGNs would have been missed by the X-ray selection 
criteria. Their conclusion is that the missing Type 2 AGNs 
are almost certainly heavily X-ray absorbed. Our sample 
shows a different behaviour. Our objects in fact do not show 
a random scatter in X-ray luminosity for a given [QUI] line 
luminosity, rather are quite aligned with the |H05| [OIII]- 
selected Type 1 sample. This can be quantified in te rms of 
the average luminosity ratios of the samples: in IhoI [oiii]- 
selected Type 2 AGNs have an average log(L x /i[oni]) of 
0.57 with a o of 1.06 dex. Our sample has instead a higher X- 



ray to [OIII] average ratio of 1.25 dex with a much smaller a 
of 0.57 dex. I nterestingly, this relation is much more similar 
to what |H05| find for the unobscured [OIII]-selected AGNs 
(log(L x /L[oiirj) = 1-59 dex, a = 0.48 dex). 

If we exclude the peculiar BL-Lac object (number 4), 
the Kendall partial correlation test shows an unexpected, 
although not strong, correlation between [OIII] and X-ray 
emissions with a significance of 99% (92% if we keep the BL- 
Lac). This correla tion is not present in the Type 2 - [OIII] 
selected sample of H05. A similar trend is instead present in 
the X-ray selected sample, where X-ray absorption removes 
all the sources in the lower right part of the X-ray - [OIII] 
plane, and in the Type 1 objects in both lH05l samples, where 
obscuration is not important. 

How can these difference be explained? The similarities 
between Seyfert 1 objects and our sample of low luminosity 
AGNs in the [OIII]-X-ray plane suggest a lack or a reduced 
incidence of molecular/dusty structures (like the dusty torus 
present in Seyfert galaxies). The effect of such a structure 
would be to absorb X-ray emission coming from the inner 
part of the AGN. This is also confirmed by the X-ray spec- 
tral analysis, which does not find any sign of absorption 
in any source of the sample (apart from number 8, that is 
also the only HEG identified in the sample) . The extinction- 
corrected [OIII] - X-ray plane, shown in Figure |4j is indeed 
similar to the uncorrected one, although due to the lack of 
high photon counts, our absorption-corrected data have big 
error bars. 

From this analysis we can conclude that, absorption be- 
ing minimal, X-ray emission could be used as a good selec- 
tion criteria for LEG sources at these low [OIII] luminosities. 
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Figure 4. The same relation shown in Figure \3\ after correction 
for absorption both in [OIII] and X-ray. Black line is the linear 
regression (object 4 excluded, slope: 0.22 ± 0.13, no evident cor- 
relation present as stated in Table [3}. Symbols are as in Figure [T] 



4.2 The radio contribution 

Our multi-wavelength dataset reveals another interesting 
fact. We find that a part of the scatter in the X-ray - [OIII] 
relation can be directly linked to the radio spectra of the 
sources. Splitting the sample into radio steep- (a > 0.3) and 
flat-spectrum (a < 0.3) sources, we find that the two sub- 
samples are still correlated in the [OIII] - X-ray plane (at a 
significance level >99% for flat-spectrum sources and of 95% 
for steep-spectrum sources, see Table [3} and their linear fits 
slopes are compatible (flat: 0.91 ± 0.09 - steep: 1.2 ± 0.5). 
But the steep-spectrum sources are systematically (~ 1 or- 
der of magnitude) more [OIII] luminous than flat-spectrum 
sources for a given X-ray luminosity. Radio emission spec- 
trum thus plays a role in the location of the source in the 
£[om] ~~ Lx plane, both for the dust-uncorrected (Figure [3]) 
and for the dust-corrected one (Figure [3]). 

One possibility is that our flat-spectrum sources occur 
in a gas-poorer environment, consequently their jets undergo 
a small interaction with the medium, so that, in the radio 
band, the very core of the AGN is the brightest part of the 
object. The flat spectrum is then explained as superimpo- 
sition of many sel f-absorbed synchrotron sp ectra generated 
by the AGN core (jBlandford fc Kon igl 1979). These sources 
will also have a fainter [OIII] emission because there is less 
gas to ionize. This picture is reinforced by the fact that 
sources with steep spectrum show on average more dust ex- 
tinction. In these sources [OIII] flux is corrected for absorp- 
tion (with the Balmer decrement method described in Sec- 
tion l3.3jl by a factor of 17±7 on average, while flat-spectrum 
sources have an average correction of a factor 2.9 ± 0.8. 

In the radio - X-ray plane (Figure [5]) there is no evident 
correlation between the two variables and the separation be- 



Figure 5. X-ray - radio (4 GHz) correlati on. Green line is from 
a sample of Low Luminosity Radio Galaxy jPanessa et aljfeoOTf) . 
Symbols are as in Figure [T] 

tween steep- and flat-spectrum sources is not present any 
longer. However, sources of our sample classified as LEG lie 
in the same location (with the im portant exception of source 
7) found bv lPanessa et all l|2007h for a sample of Low Lumi- 
nosity Radio Galaxy (green line in Figure Q extracted from 
the F R I galaxies of the 3C catalogue |Balmaverde et al] 
2006). FR I galaxies are indeed almost always classified as 
LEG, and Seyfert galaxies are usually 2-3 order of magni- 
tude more powerful in X-ray for a given radio emission. 

4.3 The fundamental plane of active black holes 

Almost all sources of our sample lie above the fundamen- 
tal pl ane of black hole activity des cribed in iKording et al.l 
(2006) (Figure left panel) and in |M03| (Figure El central 
panel). These two planes were derived using AGNs selected 
with rather different criteria. For the first using only AGNs 
where X-ray emission is believed to be jet-dominated while 
for the second explicitly excluding jet-dominated systems. 
An explanation of such behaviour in terms of underesti- 
mating the X-ray extinction from the spectra and, there- 
fore, a shift towards the left side in Figure [6] is unlikely. 
Even taking into account all difficulties in estimating the 
unabsorbed X-ray emission, the vast majority of our sources 
(14 out of 17) would need to be so heavily obscured as to 
bring their X-ray emission down by a factor ~ 100. Such a 
large fraction of heavily obscured objects would be very un- 
usual, with current estimates of the incidence of Compton 
thick AGN at these redshifts (albeit these are for somewhat 
higher luminosities) ranging between 20% and 50% (see e.g. 
iGoulding et"afll201Ch . 

Also, the possibility of overestimating radio emission 
due to boosting effect appears improbable. A boosting ef- 
fect should affect mainly sources with a flat spectrum, that 
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Figur e 6. Fund amen tal plane of acti ve black hol e relat ions. Green dashed lines shows the expected relation as found bv lKording et aD 
fcOOef) (left), bv lM03l (centre) and by I Yuan etafl J2009t) (right). Symbols are as in Figure [l] 



is likely to be connected with emission from AGNs core and 
jets, and would have affected both radio and X-ray emis- 
sions, while what we find is an excess of radio and a lack of 
X-ray emission. Furthermore, this would have brought forth 
a much larger scatter in the fundamental plane relation and 
the presence of some broadening in emission lines that are 
missing. 

In fact, LEG sources in our sample describe a tighter 
relation (again with the exception of s ource 7) than what 
expected from the fundamental plane of iMOSt This relation 
has an intrinsic dispersion of a± = 0.60, while our LEG sam- 
ple has a dispersion that is a factor 1.7 smaller (a± = 0.35, 
that goes down to a± — 0.26 without sour ce 7). The ex- 
pected dispersion of the fundamental plane of lKording et al.l 
( 2006) (a± = 0.28) is instead compatible with what we find. 
It must also be noted that part of the dispersion in our data 
certainly arises from the errors associated to (at least) our 
very weak X-ray luminosities and consequent difficulties in 
estimating absorptions. 

We find that neither lMoj nor lKording et alj jjooj ) fun- 
damental planes are able to correctly recover the relation be- 
tween radio and X-ray emissions for our sample of LLAGN. 
This is indeed expected in the first case, where the major- 
ity of jet-d ominated objec t s wer e excluded from the deriva- 
tion. In the lKording et all (|2006l ) version of the fundamental 
plane relationship instead only jet-dominated system (low- 
hard state XRB, LLAGN and FRI radio galaxies) were ex- 
ploited for derivation, nevertheless our data do not follow 
that relation either. 

Similar results were found by lHardcastle et alj (|2009h , 
they tested the fundamental plane with a set of powerful 
radio galaxies from the 3CRR catalogue whose jet and ac- 
cretion X-ray components can be disentangled using deep 
XMM spectroscopy. They found that, when considering only 
the accretion-related X-ray component, their sources lie in 
the predicted position on the fundamental plane of M03. 
Using the jet related component instead, the sources move 
to the same location where our sample of low luminosity 
objects lie. This underline the care that needs to be taken 
before using fundamental plane relations with objects that 
could be jet-domina t ed. 

In lYuan fe Cuil l|2005h an interpretation of the funda- 
mental plane in terms of ADAF-jet model is given together 
with a prediction that the X-ray emission should originate 



from jets rather than ADAFs when the X-ray luminosity 
in the 2-10 ke V band is lower th an a critical value Lx.crit 
(equation 7 in I Yuan fc Cuill2005l ). In such a situation, the 
slope of the relation between radio and X-ray luminosity is 
steeper and the normali zation higher than that predicted by 
iMOSl . lYuan et al.l l|2009l ) tested these ideas with a sample of 
22 low luminosity AGNs and found a relation very close to 
that defined by our sample (left panel in Figure [SJ) and close 
to the prediction of I Yuan fc Cuil {2005) for X-ray emission 
produced at the base of the jet. However, th e critical X-ray 
luminosity predicted by I Yuan fc Cuil ([2005), for the object 
in our sample, should correspond to log Lx.crit ~ 39.6 erg/s, 
while all our objects have higher X-ra y luminositie s. Nev - 
ertheless, it should be noted that the lYuan fc Cul (|2005l ) 
criterion is to be intended in a statistical sense and not as 
a sharp division between two populations, and the exact 
value of the critical luminosity may depend on the specific 
incarnation of the ADAF (viscosity, electron to ion heating 
efficiencies, etc.). Even so, the location of our sources in the 
radio luminosity - X-ray luminosity - SMBH mass space en- 
hance the interpretation that X-ray emission in our sample 
(at least for LEG sources) is probably related to the jet syn- 
chrotron emission rather than the Comptonization emission 
from the (radiatively inefficient) accretion flow. 



5 DISCUSSION 

As we argued in Section 14.11 at least for LEG sources, the 
absence or a reduced incidence of dusty /molecular gas may 
explain the [OIII] - X-ray correlation and the unusual miss- 
ing signs of X-ray absorption in all LEG X-ray spectra (signs 
of absorption are only found in the HEG source). Moreover, 
part of the scatter in the observed [OIII]-X-ray relations is 
directly linked to the radio spectra of the sources, and we 
have argued in Section 14.21 that our flat-spectrum sources 
live in a gas poorer environment, and consequently their 
jets undergo a small interaction with the ambient medium. 

Independent evidence for an intrinsic difference in the 
amount of cold/obscuring ma terial in low l u mino sity AGN 
has been gathered recently. Elit zur fc Hoi |2009l ) showed 
that, in a nearly complete sample of nearby AGN, the broad 
line region disappears at low intrinsic luminosities, and inter- 
preted this fact within the disk-wind scenario for the broad 
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line r egion (BLR) and toroidal obscuration in AGN (|Elitzurl 
2008). If the BLR constitutes the inner, du st-free part o f 
the cold molecular obscuring torus (see e.g. iNetze r 2008), 
then a similar behaviour should be expected al so for other 
indica tors of intrinsic nuclear absorption. Indeed pTr an et all 
(2011), using optical spectropolarimetry on a set of 3 low 
luminosity Type 2 AGNs, found no si gn of broad l i ne re - 
gion nor obscuring torus. Furthermore, iBurlon et al.l (|2010T l 
found for the first time tantalizing evidence that the fraction 
of obscured AGN in a complete, hard- X-ray selected sample 
of nearby AGN, declines at Lx < 10 42 erg/s. 

In Section T4. 31 we pointed out a diffe rent b ehaviour be- 
tween LEG sources and those defining the lM03l fundamental 
plane, and this is explainable with a different X-ray produc- 
tion mechanism between the population used to construct 
the fundamental plane and our sample. In the fundamen- 
tal plane, X-ray emission is assumed to be generated by a 
(radiatively inefficient) accretion flow while, in LEG objects 
presented here, a much weaker (about 2 order of magni- 
tude) X-ray emission could be generated by the jet itself. 
The tight correlation (a± — 0.35) between radio and X-ray 
emission for LEG objects enhan ces this conclusion , and is 
in agreement with the analysis of I Yuan et all (|2009l L Future 
studies on the relation between radio and X-ray emissions 
and the object mass in a bigger sample can unveil possible 
new relations similar to the fundamental plane but tuned for 
low luminosity, jet-dominated objects. Based on the results 
presented here, s uch a plane will lie above the fundamen- 
tal plane found in lM03l and will have a smaller scatter as a 
consequence of the tight relation between radio and X-ray 
production in such kind of sources. 

With respect to the analysis performed in iMOl we re- 
mark here that the heterogeneity of their sample selection 
prevented any uniform classification of the X-ray spectra in 
terms of the amount of obscuration present. Based on the 
work presented here, it is likely that a large part of the scat- 
ter in the original fundamental plane relation could be due 
to the inclusion of both gas poor and gas rich sources in the 
original analysis. 

Taken together, these considerations open the possibil- 
ity to use X-ray emission as an unbiased estimator of source 
total power for low luminosity AGN, provided robust in- 
dicators of a lack of cold gas in the observed systems is 
available. [OIII] emission, although useful in the detecting 
and classification procedure, can be biased by the amount 
of dust present in the AGN surroundings, as we can infer 
by different position of steep- and flat-spectrum sources in 
Figure [4] 



6 CONCLUSIONS 

In this paper we analysed a set of sources extracted from the 
SDSS catalogue. For these sources we have measured Chan- 
dra X-ray flux, radio luminosity in two bands from the VLA 
and optical spectral information extracted from the SDSS 
database. Our sources span a wide range in optical [OIII] 
and radio luminosities but they have very similar SMBH 
masses (~ 10 s M Q ) and are all located over a narrow red- 
shift range (0.05 < z < 0.11). 

Although our sample spans a wide range of [OIII] and 
radio luminosities, we found that 15 out of 17 sources are 



classified as LEG. The sample shows a correlation on the 
uncorrected X-ray - [OIII] luminosity plane that is not seen 
in sa mples of Type 2 AGN with higher [OIII] luminosity 
(H05). An absent or reduced obscuring torus can explain 
this relations and different amount of gas in the AGN sur- 
roundings can account, to first order, for its scatter. This 
last property is well connected to the radio spectral index. 
This has important implications for selecting criteria in up- 
coming low-luminosity AGN surveys. For LEG samples, the 
lack of X-ray obscuration means that X-ray emission can be 
used as a good selection quantity with no obscuration se- 
lection biases. Nevertheless, radio core emiss ion is stil l the 
easiest way to identify this class of objects (|Hol I2008T ) . al- 
though more information (i.e. from optical spectra) is still 
required for object classification. 
Summarizing, our conclusions: 

• At low luminosities (Lx < 10 42 erg/s) a selection cri- 
teria based on [OIII] emission is not required in order to 
minimize the loss of of obscured sources. At low luminosity, 
X-ray and [OIII] emissions are correlated and there are no 
evident signs of X-ray absorption. 

• We found that all steep spectrum sources have around 
an order of magnitude less [OIII] emission than flat- 
spectrum sources with similar X-ray luminosity. We argue 
that the amount of interstellar medium available for inter- 
action in the proximity of the AGN is responsible for this 
distinction in the radio spectrum and for the difference in 
the production of the [OIII] line. 

• In the radio - X-ray relation (and so in the fundamental 
plane of active black holes) , LEG sources have a smaller dis- 
persion and a reduced X-ray production compared to other 
types of AGN (i.e. Seyfert galaxies and quasars). This can 
be explained associating X-ray and radio production to the 
same physical mechanism, i.e. synchrotron emission from the 
base of the jet, i n very good agreement with recent inves- 
tigation (see e.g.. lHardcastle et al.ll2009T ) although obtained 
from sample selected with different criteria. 
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